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ABSTRACT 

We use the mean-field dynamo equations to show that an incoherent alpha 
effect in mirror-symmetric turbulence in a shearing flow can generate a large 
scale, coherent magnetic field. We illustrate this effect with simulations of a 
few simple systems. In accretion disks, this process can lead to axisymmetric 
magnetic domains whose radial and vertical dimensions will be comparable to 
the disk height. This process may be responsible for observations of dynamo 
activity seen in simulations of dynamo-generated turbulence involving, for 
example, the Balbus-Hawley instability. In this case the magnetic field strength 
will saturate at ~ (h/r)'^ times the ambient pressure in real accretion disks. 
The resultant dimensionless viscosity will be of the same order. In numerical 
simulations the azimuthal extent of the simulated annulus should be substituted 
for r. We compare the predictions of this model to numerical simulations 
previously reported by Brandenburg et al. (1995). In a radiation pressure 
dominated environment this estimate for viscosity should be reduced by a factor 
of {Pgas / PradiationT ^ue to magnetic buoyancy. 



1. Introduction 



Understanding the transport of angular momentum in accretion disks is one of the 
basic challenges in modern astrophysics. The traditional approach (|Shakura, fc Sunyaev 



1973|) is to assume that accretion disks are characterized by an effective viscosity, arising 



from an unspecified collective process, given by assCsh, where Cg is the local sound speed, 
h is the disk half-thickness, and ass is a constant of order unity. More recently, there has 
been the realization ( [Balbus fc Hawley 199 1| ) that a previously discovered magnetic field 
instability in a shearing flow ( |Velikhov 1959| , phandrasekhar 1961| ) will act to produce a 
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positive angular momentum flux in an accretion disk. This has given rise to two separate, 
but related claims. The first is the proposal that this is the dominant mechanism of 
angular momentum transport in ionized accretion disks. The second is the proposal that 
this instability, by itself, leads to a turbulent dynamo which drives the magnetic field into 
equipartition with the ambient pressure, i.e. Va ~ c^, where Va is the Alfven speed in the 
disk. The growth rate for this dynamo is usually taken, following the original claim of 
Balbus and Hawley, to be ~ Q. Since the dimensionless 'viscosity', ass, is ~ (Va/cs)^, this 
proposal requires that ass be a number 'of order unity'. However, we need to interpret 
this phrase generously. In numerical simulations (e.g. [Brandenburg, Nordlund, Stein, & 



Torkelsson 1995| ) ass is less than 10 ^ both because the magnetic pressure saturates at a 



fraction of the gas pressure, and because the off-diagonal components of (BB) are a fraction 
of Bl 

Three dimensional simulations of the Balbus-Hawley instability have been performed 
by a variety of researchers, with and without imposed vertical magnetic flux, and with and 
without vertical disk structure ( [Hawley, Gammie, fc Balbus 1995"a| , Prandenburg, Nordlund 



Stein, fc Torkelsson 1995 , Stone, Hawley, Gammie, fc Balbus 1995[ , and Hawley, Gammie" 



fc Balbus 1995b[) . We note in particular Brandenburg et al. (1995) in which no net flux 
was imposed on the computational box, and vertical disk structure was included. In this 
simulation, as in the others, there was an initial rise in the magnetic energy density at 
a rate ~ Q. At the end of this phase the system had not yet lost memory of its initial 
conditions, but after a somewhat longer time, which may be as long as a few dozen rotation 
periods, the simulation asymptotically approached a final state with Va ~ c^. The approach 
to this state was characterized by the appearance of a large scale field which underwent 
spontaneous reversals at irregular intervals of tens of rotational periods. 

Interestingly, the presence of a large scale coherent field does not seem to be due to 
an a — Q dynamo, because the relative helicity is just a few percent. Conventional a — Q 
dynamo models rely on a nonzero aee component in the helicity tensor (not to be confused 
with the dimensionless viscosity, written here as ass) to produce a large scale coherent 
field. The presence of an initial rapid rise is less surprising, since imposing a uniform large 
scale magnetic field in a turbulent medium results in the formation of intermittent magnetic 
field structures and a consequent rise in the magnetic energy density at the turbulent eddy 
turn-over rate. In addition, there is evidence ( [Meneguzzi, Frisch, fc Pouquet 1981[ ) that 



turbulence in a conducting fluid can generate a modest and highly disordered magnetic 
field even in the absence of an imposed global field. Both of these effects are probably due 
to the ability of of symmetric turbulence to produce a negative effective diffusion coefficient 
( [Mottatt 1978[ ) and they both can be relied open to contribute to the growth of the high 



wavenumber component of the magnetic field. On the other hand, the slower relaxation rate 



- 3 - 



seen after the initial rise is correlated with changes in the large scale field and is presumably 
an indicator of the large scale dynamo growth rate. Since the turbulence is sustained by an 
instability of the large scale field, its ability to generate such a field is critically important. 

The saturation level of the magnetic field in these simulations also leads to some 
puzzling questions. The claim that the Balbus-Hawley instability saturates when the 
magnetic pressure is comparable to the ambient thermal pressure, and that the dimensionless 
viscosity has some approximately fixed value, is difficult to reconcile with attempts to model 
accretion disks in compact binary systems. Successful models of dwarf novae outbursts and 
X-ray transients ( |Smak 1984a| , |Smak 1984b| , [Meyer fc Meyer-Hofmeister 1984| , [Huang & 



Wheeler 1989| , pVlineshige fc Wheeler 1955| , and, more recently, [(Jannizzo 19941 ), as well as 



the distribution of light in quiescent dwarf novae disks ( jMineshige, fc Wood 1989D all imply 



that the dimensionless viscosity, ass, varies spatially and with time. These variations are 
consistent with ass oc (h/r)"', where n is a constant lying somewhere between 1 and 2. 
Recent work ( Pannizzo, Chen, fc Livio 1995|) on X-ray transients suggests that n may be 
close to 1.5. Here we note only that any value of n appreciably different from zero confiicts 
with claims for a universal value of ass- 

This difficulty can be resolved in several different ways. For example, we might claim 
that magnetic instabilities dominate ass only at low temperatures and that some other 
process, e.g. convection, dominates at higher temperatures. This idea faces two major 
objections. First, it explains only some of the phenomenological evidence favoring a varying 
ass- Second, attempts to model the vertical structure of dwarf novae disks invariably 
conclude that such disks are convective during quiescence, when ass is small and stably 
stratified during outburst, when ass is relatively large (for a recent discussion of the 
conditions necessary for convection in partially ionized accretion disks see Cannizzo 1992). 
This implies that convection could explain the rise in ass only if it acts to suppress angular 
momentum transport, rather than enhance it. Alternatively, one could appeal to the 
temperature dependence of the resistivity to account for this effect, although the effective 
resistivity of the simulations is, in any case, many orders of magnitude larger than in real 
disks. A more promising notion is that one might ascribe the rise in ass to the greater 
thermal conductivity of disk in the hot state, although the rationale for this is not yet clear. 
Finally, one might simply conclude that all the phenomenological models are wrong, for a 
variety of reasons, a viewpoint which is difficult to dismiss given the large uncertainties 
faced in modeling accretion disks. 

In this paper we will explore a new turbulent disk dynamo in which the turbulence is 
not assumed to lack mirror symmetry in the vertical direction. The dynamo effect arises 
from the fact that in a system of finite size the mean square helicity, and the instantaneous 
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spatially averaged helicity, is still nonzero. We will see that this leads to the existence of a 
modified a — Q dynamo, in which the large scale organization of the magnetic field comes 
from the existence of a large scale shear. Of course, a real disk has vertical structure, which 
breaks the vertical symmetry and allows for the possibility of a nonzero average helicity. 
However, we will show that the incoherent dynamo mechanism will be particularly effective 
in simulations of limited azimuthal extent. Indeed the simulations of Brandenburg et al. 
(1995) show that while the relative helicity is small (less than a few percent), there is still a 
dynamo effect leading to the generation of large scale fields. Moreover, the saturation level 
of the magnetic field, and the consequent value of ass, turn out to depend sensitively on 
the ratio h/r. If this dynamo is the only effect to arise from the turbulence induced by the 
Balbus-Hawley effect, then it is relatively simple to reconcile the dynamo activity seen in 
simulations with phenomenological models of accretion disks in compact binary systems. 

In §2 we discuss the conceptual basis of an incoherent dynamo in a turbulent shearing 
medium and estimate the growth rate. In §3 we apply this to accretion disks and show 
that the incoherent dynamo gives a positive growth rate only for axisymmetric magnetic 
domains. We estimate the saturated state of the field and discuss our results in light of 
numerical simulations of magnetic fields in a Keplerian disk. In §4 we summarize our results 
and their implications for astrophysical disks and numerical simulations of such disks. 



2. The Incoherent Dynamo 

In a highly conducting medium the magnetic field obeys the induction equation 

dtB = V x{Vx B), (1) 

where we have neglected ohmic diffusion. Ultimately this term is important in allowing 
reconnection and smoothing. Here we assume that these processes take place at a rate 
determined by turbulent processes. The usual approach to dynamo theory is to define the 
response of the large scale magnetic field to small scale motions as b and to derive its effects 
on the large scale field by substituting b back into the right hand side of equation (p. For 
an incompressible fluid this yields 

dtBi = eijkdj{akiBi) + djDjkdkBi - djDikdkBj, (2) 

where eijk is the Levi-Civita tensor and 

aj,i = {e,,,VA fv,{t')dt'), (3) 
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and 



D,k = (V,- j Vk{t')dt') 



(4) 



The first term comes from tlie stretcliing of large scale field lines by the local turbulence. The 
tensor a^i describes the twisting of large scale field lines into a spiral shape. Reconnection 
between adjacent spirals produces a large scale field component at right angles to the 
original field line provided that either the degree of twisting or the large scale magnetic field 
strength varies in the third direction. The second term is the usual diffusion term, modified 
by the presence of the third term. 

We can see from equation (Q) that each component of aij has either a factor of or 
dz- If the local velocity field is mirror symmetric, in the sense that its statistical properties 
are unchanged under the transformation z —z, then the time and space averaged value 
of aij vanishes. This poses a significant but, as we will see, not insurmountable, obstacle to 
a successful dynamo. 

Another problem is that equations (0), (H) and (^) are usually defined kinematically, i.e. 
the velocity field is assumed to be imposed on the magnetic field. Once the magnetic field 
becomes sufficiently powerful it will modify the flow, which is usually taken into account 
by including a correction term proportional to B"^. However, in a Keplerian shearing flow 
the magnetic field will be unstable and the resulting turbulence will be directly correlated 
with the magnetic field. Nevertheless, as long as we define V in terms of the motion of the 
magnetic field lines equation (j^) will remain valid, if difficult to solve. Here we will define 
our results in terms of the properties of aij and Dij regardless of their ultimate source. 

In a Keplerian disk the dynamo equations can be simplified as 



where Vt oc r~^/^ is the rotation frequency, Vh is the buoyant velocity of the magnetic 
field lines relative to the surrounding fluid, and Br and Bq are the radial and azimuthal 
components of the magnetic field. Equations (§) and differ from equation (H) in that 
we have allowed for the presence of global shearing, and magnetic field line buoyancy. In 
addition, we have assumed that the diffusion matrix is diagonal, and dropped the effects 
of helicity on the evolution of Sg, given that the shearing of Br should dominate such 
effects. Also, we have retained only the age term in equation (j^) since the critical feedback 
term in the dynamo equations involves generating radial magnetic flux from the azimuthal 
component of the field. Finally, given that we are interested in applying these equations to 




(5) 



and 
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accretion disks whose thickness is a small fraction of their radius, we have assumed that 
vertical gradients will dominate over radial gradients. 

Now let's assume that the turbulence is symmetric under z — >• — z so that (aee) = 0. 
Although this eliminates any coherent helicity, the value of {B^) can still increase in a 
random walk. Ignoring diffusion and buoyancy we see that the formal solution for is 

Br = f -d,{aee{t')Be{t'))dt'. (7) 

Now by hypothesis, agg is uncorrelated over time scales greater than some eddy correlation 
time Teddy If the radial magnetic field is undergoing a random walk, then it will usually 
be far enough away from zero that it will not change sign every eddy correlation time. 
Since B^. drives Bg through coherent shearing, this implies that the correlation time for 
Br and Bg is much greater than Teddy Consequently, we can consider the integrand in 
equation as consisting of a rapidly varying factor, agg, multiplying a slowly varying 
function. Multiplying equation (0) times equation (|) and ignoring diffusion and buoyancy, 
as before, we see that the integral in equation (|^ is correlated with agg only over the last 
eddy correlation time Teddy Consequently, we can replace the integral in the product with 
-dz{agg{t)Bg{t))Teddy This implies 

d,{Bl)^Kl^flTeddy{Bl), (8) 

where is the number of independent turbulent eddies in a magnetic domain, Kz is 
the vertical wavenumber of the magnetic domain, and {a^g) is the mean square helicity 
associated with a single eddy. In general this will be of order where Vt is the root mean 
square turbulent velocity. Since B^. is being driven incoherently we can expect it to undergo 
frequent reversals. In between such reversals the shearing of the field will drive Bq sharply 
upward. From equations and (^) we see that the correlation time of the radial magnetic 
field and the growth time of the magnetic field are comparable and given by 



corr growth 

We note that Teor-r has to be greater than Teddy in order for this estimate to be internally 
self-consistent, i.e. the magnetic field must be correlated over longer times than the 
turbulence itself. Since a field reversal in Bg requires that Br not only reverse its sign, but 
maintain it long enough to push Bg through zero, it is clear that the correlation time for Bg 
may be somewhat larger than the correlation time for Br- We will return to this point later. 

By itself this argument does not show that a succession of random twists in a shearing 
background can drive an exponential increase in the magnetic field. We need to show that 



K^zVSn^T\ 
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the growth experienced between field reversals dominates over the abrupt cancellation of 
the field as Br reverses itself. We also need to show that our estimate of the growth rate 
given in equation (j^) will dominate over turbulent diffusion for some range of magnetic 
domain sizes. 

We can test the assertion that a series of random changes in B^. can drive a dynamo by 
constructing a simple toy model of the process, which ignores the spatial structure of the 
field, but includes its dynamical evolution. Assuming that age has a stochastic component 
and ignoring buoyancy we can rewrite equations (|) and (^) as 

dtBr = ivit) - oicoh)Be - DBr, (10) 

and 

dtBe = --VlBr - DBe, (11) 

where rj{t) is a stochastic variable with a correlation time Teddy and acoh is the coherent 
component of dzCtee- Here we have subsumed spatial derivatives into the definitions of r] and 
D and ignored the —aeedzBe term which would normally appear in the mean-field dynamo 
equations. We have also assumed that turbulent damping is the same for each component 
of the magnetic field, which is not generally true, but simplifies the analysis without losing 
any essential physics. Equations ([lO|) and (11) can be rewritten in a more convenient form 



by defining A = {Br/ Be). Then 

dtA = r]{t) - aeok + ^nA^, (12) 

and 

dtlnB^ = -3nA-2D. (13) 

The magnetic field will grow exponentially if (A) is negative and —3Q{A) > 2D. 

We can find (A) by solving equation (|12D in terms of an unnormalized probability 
distribution function P{A) and evaluating 

^ IZoAPiA)dA 
^ ^ - IZoPiA)dA- ^'^^ 

The distribution function P{A) satisfies the equation 

dA{AP{A) - {r]^)TeddydAP{A)) = 0, (15) 



or 



P{A) (^^QA^ - aeoh) - {v')reddydAP{A) = {v')reddy, (16) 
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where we have taken advantage of the unnormahzed nature of P{A) to set the constant of 
integration to {rf)Teddy Equation ( ]T6| ) can be solved to yield 



P{A) = exp 



nA^ - 2acohA 



Consequently, 



'2{7]'^Teddy) 



2{v'nd,y) \ I-oo dy ly^ dsy expb^ - - ^{y - s)] 



exp 



2{7f Teddy) 



dr. 



(17) 



{A) 



I-oo dy ds exp[?/3 - s3 - 7(?/ - s)] 

1/3 



where 

{rj^Teddy) \ VL ) 

Equation (0) can be rewritten by defining new variables w = y + s and x = s — y and 
integrating over w. We obtain 



^ -1 ( 2{r^'reddy) \ x^l^ eM<l - i)]dx 

^ ^ 2 1^ ) /-a;-V2exp[x(7-f)]rfa;- 

When 7 is small we can expand e"''^ ^ 1 + 7W and obtain 

^2 



(20) 



which implies that 



. iV^reddy) J 



1/3 



2D. 



(21) 



(22) 



In other words, the magnetic field will grow exponentially roughly as fast as the estimate 
given in equation (H). This will be suppressed by turbulent diffusion only when the damping 
rate due to diffusion is comparable to the growth rate. 

The existence of an incoherent dynamo emerges from the fact that the distribution 
function P{A) given in equation (^) is biased towards negative values of A. This bias 
comes, paradoxically enough, from the coherent, positive definite term in equation ([12|) . 
When A is sufficiently positive it evolves deterministically through +00 into negative 
values. (Actually, doesn't change during this phase. This deterministic trajectory is 
merely a field reversal for Bg.) The end result is that whenever A becomes large and 
positive it rapidly switches to being large and negative. Ultimately, the sign of the bias is 
determined by the sign of dr^. The frequency of such field reversals is given by examining 
the probability distribution at large A when the evolution of P{A) is deterministic. If we 
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define t{A) as the time it takes for tlie field to move from some large positive value of A to 
A = oo then from equation (|12|) we see that 

r{A)-' = ^nA, (23) 

where we have neglected acoh since for A sufficiently large its effects can be ignored. The 
field reversal rate is just the limit of this rate times the statistical weight of the distribution 
between A and oo. In other words, 

-.-i= limr(A)-ir^^. (24) 

Substituting equations (|1^) and (^) into this result, and making the change of variables, 
as before, to x and w we have 

r-l = hm 3^^ /2A(o/2(,^.....))V3/o-exp ^ , ^^^ .^ 

^--2 JZoirexp[w{-f-\{w^ + 3x^))]dwdx 

Both the numerator and the denominator can be simplified by integrating over x to obtain 

'''' U/ /o°°w;-i/2exp[u;(7-«;V4)]ciw' ^ ' 

When 7 is small this becomes 

r-i = 0.53((r7V,rfrf,)^]2)l/3(l _ 0.537), (27) 

i.e. a rate which is roughly half the e-folding rate for the magnetic field energy. 



When acoh is large and positive we can evaluate the integrals in equation ([20| ) by 
expanding around the maximum of x(7 — x^/r). We obtain 

(^)--(^) , (28) 

SO that 

dt\nBl^2{^acoh^ -2D, (29) 

which is the expected result for a coherent a — Vt dynamo. In this limit the magnetic field 
reversal rate becomes 



A.lacoh 



(30) 
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As expected, field reversals are exponentially suppressed as we go to the usual a — Q 
dynamo. Given acoh > then as acoh becomes significant we expect it to enhance the 
dynamo growth rate and reduce the rate of spontaneous field reversals. 



When acoh is large and negative we can evaluate equation (^0]) by integrating the 
denominator by parts and remembering that the bulk of the contribution to the integral 
comes from w < —l/'j so that w'^ <^ —7. We obtain 

(A) ^ i^!£^, (31) 



and 



In this limit field reversals occur at a rate given by 

= 0.78i-acohnY^'. (33) 

We note that in this case the coherent component of the helicity does not completely 
shut off the incoherent dynamo, even though by itself it is incapable of driving a dynamo. 
Instead we find that as I7I increases past one the dynamo growth rate decreases inversely 
with I7I. Eventually, turbulent diffusion will suppress the dynamo. In the limit where 7 is 
of order —1 we anticipate that the dynamo growth rate will be less than expected from the 
incoherent dynamo alone and the rate of field reversals will be larger. 

These analytic results have the advantage of being based on a solvable model, but 
do not include the effects of spatial structure or saturation. It is therefore instructive 
to consider the combined effects of a random electromotive force and shear in a one 
dimensional model. We consider the mean field equations for a uniform disk with Keplerian 
rotation and half-thickness H, 

dtB, = -d,{aBg) + DtdlB,, (34) 

dtBe = -^ilBr + Dtd^Be, (35) 

with —H<z<H and Br = Bg = at z = ±H. We first consider the incoherent a-effect, 
so we take a to be random in space and time. When the rms value of a is large enough 
we find self-excited solutions that grow without bound. In reality there must be some 
quenching mechanism, which we model using 
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where AA is a random function in space and time with zero mean and an rms value of unity. 
Without loss of generality we put H = Q = Dt = I. 

In fig. 1 we plot contours of the Bg field in a space-time diagram for a dynamo 
number aoQH^ / of 10^. (The critical dynamo number for dynamo action depends on 
the coherence time and length scales, A and r, respectively. In the present case we adopt 
A = 0.05 and r = 0.002 and find the critical dynamo number to be around 2000. At this 
dynamo number the ratio of the growth rate given in equation to Df/H"^ is ~ 7.) The 
remarkable result is that the Bg field shows a great deal of spatio-temporal coherence 
with variations comparable to the diffusion time and diffusion length. Experiments with 
different dynamo numbers suggest that the degree of coherence is more pronounced for 
larger dynamo numbers. 

In order to isolate the effect of a spatially incoherent a-effect we now investigate a 
model with a steady a-effect of the form 

a = ao sin(n7r2;). (37) 

For large values of n, the critical dynamo number is proportional to n"^. Thus, although 
the rms value of the a-effect is unchanged, the dynamo becomes harder to excite if a is 
chopped into many domains of different sign. The magnetic field is steady, and the radial 
component is of alternating sign. However, more surprisingly, the toroidal magnetic field 
has the same sign for all values of z, see fig. 2. This is very similar to the simulation of a 
random incoherent a-effect mentioned before. There is one difference in that the magnetic 
field shows global reversals in time when the a-effect is incoherent in time. 

Finally, we note that in order for the magnetic field to grow the growth rate given in 
equation (P) has to be greater than the dissipation rate. In general the dissipation rate will 
depend on the wavenumber of the magnetic domain as K"^, while the growth rate goes as 
(K"^ /NY^^. Clearly whether or not there is a self-excited dynamo will depend in large part 
on the geometry of the fluid. 



3. The Incoherent Dynamo in Accretion Disks 

A Keplerian accretion disk with a root mean square Alfven speed of Va will be subject 
a local instability first described by Velikhov (1959). Its pivotal role in transporting angular 
momentum outward in accretion disks was recognized later ( [Balbus fc Hawley 1991| ). In 
the context of accretion disks this instability is normally referred to as the Balbus-Hawley 
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instability. Its maximum growth rate is of order f2, and occurs at an azimuthal wavelength 
of ~ Va/^. In three dimensions the instability saturates in turbulence with a typical 
turbulent velocity comparable to Va and a typical eddy size of ~ Va/^- This turbulence 
is not expected to be isotropic, but the typical eddies are expected to have axis ratios of 
order unity, which in this context means only that no axis should be more than an order 
of magnitude larger than another (|Vishniac, fc Diamond 1992|) . Numerical simulations 



( Prandenburg, Nordlund, Stein, fc Torkelsson 1995|) indicate that the the azimuthal scale of 
the typical eddies is several times the vertical and radial scales, which is expected in light 
of the large local shear. The azimuthal velocity is also larger, although only by a factor 
of roughly two. Neglecting such factors, these scaling laws imply a turbulent diffusivity of 
~ Vj^/Q. The turbulence largely suppresses the Parker instability and the typical buoyant 
velocity of the magnetic field is of order V^/cg ( |Vishniac, fc Diamond 1992| , [Vishniac^ 
|1995b| ), where Cg is the local sound speed. The angular momentum flux induced by the 
turbulence is approximately (VeVr) ~ V| which implies a dimensionless viscosity ass of 
order (Va/cs)^. Since hQ ~ Cg this implies that magnetic flux is lost from the disk at a rate 
which is some fraction of order unity times ass^- 

It is by no means obvious that in real disks this turbulence possesses the kind of 
symmetry that would make (age) = 0. On the other hand, calculations done without vertical 
structure or any imposed large scale field (Hawley, Gammie, fc Balbus 1995E) give results 



which are qualitatively similar to calculations which include vertical structure ([Brandenburg 



Nordlund, Stein, fc Torkelsson 1995|) . By construction the former calculations must be 



symmetric under the transformation z —z even though the latter are not. We can in 
principle estimate agg using data from the simulation of Brandenburg et al. (1995). From 
equation (|^) it is clear that a time integration has to be carried out. However, video 
animations of those data suggest that the life time of turbulent eddies is shorter than the 
life time of magnetic structures which, in turn, is shorter than the eddy turn over time. 
In other words, the Strouhal number (e.g. Krause & Radler 1980) is small. As a rough 
approximation we may therefore replace the time integration by a multiplication with a 
relevant time scale. We adopt the natural time scale which is sufficient since we are 
only interested in relative variations. We adopt volume averages and note that because of 
the periodic boundary conditions in the toroidal direction, (KV^.e) = — (V^K.e), so we can 
compute 

aee^-{VrV,,e)^-\ (38) 
r 

In fig. 3 we plot the evolution of agg using the data from run C of Brandenburg et al. 
(1995), which has now been carried out for an additional 200 orbits, see also Torkelsson 
et al. (1996). This average was computed for the upper half plane of the simulation. We 
note that agg is positive, in agreement with the expected effect for bubbles that expand as 
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they rise in a Keplerian disk. However, the sign of agg suggested by the correlation between 
the azimuthal magnetic and electric fields is negative ( [Brandenburg, Nordlund, Stein, ^ 



Ibrkelsson 1995| ). The source of this discrepancy is not yet clear. In any case the spatially 



averaged helicity shows large variations from its long term average, although the variations 
in the electromotive force are much larger. 

The size of the fiuctuations in the electromotive force, as well as the persistence 
of the dynamo in the absence of any z symmetry breaking, implies that any preferred 
helicity resulting from vertical structure is not strong enough to completely dominate the 
simulations. In what follows we will assume that real disks lack any significant (agg). 
At a minimum our results can be taken as demonstrating that there is an incoherent 
dynamo operating in the simulations, and in real accretion disks, whose effects need to be 
understood, and cleanly separated from any other dynamo mechanisms at work. 

Let's consider a magnetic domain characterized by the wavenumbers {Kr, Kg, Kz). 
Ignoring the anisotropies in the turbulence we find that the number of turbulent eddies per 
domain is roughly ~ {KrKgKzV^)^^Q'^ . Consequently the growth rate for the dynamo is 

'dynamo ~ I ^2 j ' W"^/ 

However, in a shearing environment we aren't free to specify and Kg separately. The 
shear implies a minimum Kr for any Kg since in a time ~ Tdynamo the shear will increase Kr 
by an amount {3/2)KgQTdynamo- If we choose a value of Kr above this minimal value then 
Tdynamo will go as K^^ while the dissipation rate scales as K^. Clearly our chances for a 
successful dynamo will be maximized by taking Kr ~ KgVLTdynamo- This gives us 



1/2 



'^dynamo ~ ^ ^3 j " ('^^) 

This analysis only makes sense in the limit where the magnetic domains encompass at least 
one eddy, or K^Va < ^ and KrVA < The dissipation rate is roughly 

dissipation \ z ' rl ^ ' \ / 



By comparing equations (^) and (|41|) we see that the incoherent dynamo is incapable 
of generating non-axisymmetric large scale magnetic fields. The dissipation rate of such 
domains exceeds the generation rate for all domain sizes greater than a single eddy. 

In a real disk the number of eddies in a magnetic domain does not increase indefinitely 
as Kg —>■ 0. The finite circumference of the disk implies that for axisjTiimetric domains 

N — (42) 

KzKrVr ^^'^ 
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Consequently, we can rewrite equation (^) as 



'dynamo I 



^dynamo ~ ( „o2 I 

(43) 



At a fixed wavenumber, and therefore at a fixed dissipation rate, tliis rate is maximized for 
Kz = KrS^^^. Assuming tliis ratio we see tliat tlie dynamo growtli rate for axisymmetric 
domains goes as K^^^, wliicli implies tliat at some sufficiently small K the dynamo will 
work. More exactly, the incoherent dynamo caused by the Balbus-Hawley instability will 
drive an increase in the magnetic field strength if 

< (44) 
tVa ^ ' 

In other words, the incoherent dynamo only works for 

Va < (45) 

Ultimately Kz is limited by the height of the disk, i.e. K^h > 1. Moreover, as we 
approach this limit the buoyant loss of magnetic flux becomes significant. The buoyant loss 
rate from a single magnetic domain goes as 

nlly,^,r^KMr-Kz^, (46) 

SO when K^h 1 buoyant losses are as important as turbulent diffusion. Of course, the 
only limit on the radial extent of a magnetic domain is Kj.r > 1, but lowering K^. past 
will lower the growth rate without affecting the dissipation rate. From equation (^) we 
see that the magnetic field associated with scales of order the disk thickness will be the 
strongest and will be given by 

Va--Cs-. (47) 
r 

This in turn implies that the dimensionless viscosity associated with this dynamo mechanism 
is ^ 

ass - (^) - (^) . (48) 

We expect this scaling law to hold only in the limit h <^r. As — > Cg corrections of order 
{Va/cs) will become important in our formula for buoyancy. Since the saturation limit for 
the magnetic field involves the small difference between the growth rate dependence on Va, 
which has an exponent of 5/3, and the buoyant loss rate dependence, which goes as K4, 
we expect the saturation strength of the magnetic field to be extremely sensitive to such 
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corrections unless <^ c^. We also note that the disk radius enters into this result only 
through its role as the circumference of an annulus. Computer simulations typically involve 
a short arc in place of a full annulus. In this case the azimuthal length of the simulation 
has to be used in place of r in equation (HSI). 



The rate of spontaneous magnetic reversals expected in the absence of any coherent 
component to age is comparable to dynamo growth rate. However, while current simulations 
seem to show a significant reversal rate in the presence of vertical structure ( [Brandenburg 



Nordlund, Stein, fc Torkelsson 1995|) , in its absence the field can evolve for 100 orbital 



times without reversing ( [Ibrkelsson, Brandenburg, Nordlund, fc Stein 1996| ). The exact 
relationship between the dynamo growth rate and the field reversal rate dependent on the 
particular model for the process, and the zero dimensional model used in this paper may 
well overestimate the rate of spontaneous field reversals. Nevertheless, such reversals are 
an intrinsic part of the model, and should occur if the simulation is run for several growth 
times. The sharp rise in the field reversal rate when vertical structure is included suggests 
that a significant coherent aee is present in the such simulations. (In the zero dimensional 
model this would argue for a negative vertical gradient in the coherent aee- The situation 
is less clear in three spatial dimensions.) In order for this helicity to allow the buildup of a 
coherent field in these simulations, as well as in accretion disks, it has to scale with the local 
rms turbulent velocity more steeply than the square of the incoherent dynamo growth rate, 
or {Vt/ CsY^^^ ■ If this helicity is due to the Parker instability and if, as has been argued 
elsewhere ([Vishniac, fc Diamond 1992|) , the Balbus-Hawley instability reduces the Parker 
instability to vertical motions of order V^/c^, then we can estimate the magnitude of the 
helicity as 

VrkeV:,nuoyanU (49) 

where Tbuoyant is the correlation time for these buoyant motions. Since shearing imposes the 
requirement that kgQ < KTbuoyant ^'^d since these motions are approximately incompressible, 
i.e. kr-Vr ~ kzVz, this gives a helicity less than 



Vz'kz/n ~ (50) 



If the coherent helicity has this dependence, then it becomes important only as the dynamo 
saturates due to turbulent mixing and buoyancy. In this case it will not suppress the 
incoherent dynamo in simulations with smaller {h/r), or in real accretion disks, but it will 
remain significant in the saturated state. 

The fact that the buoyancy does not significantly enhance the loss of magnetic flux 
is a critical element in the derivation of equation (^8]). Consequently, environments that 
increase magnetic buoyancy will saturate at much lower field strengths. As an example we 
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can consider magnetic flux tubes in a radiation pressure dominated environment. In this 
case we have ( |Vishniac 1995b| ) 



P T/2 

-'^radiation ''A 



Vb — . (51) 

^ gas 



Combining this resuh with equation (^Hf) for Kr ~ Kz ^ h ^ yields 



2 



'-^s / V ^radiation - 



The exponent given in equation (|52D is perhaps a bit large in comparison to the value 
suggested by the phenomenology of disks, but this is a considerably less serious problem 
than if it were too small. Competing dynamo mechanisms and/or hydrodynamic angular 
momentum transport mechanisms could be driving ass up. How does this compare to 
other sources of viscosity in disks? The result given in equation (|52D has an extremely 
uncertain coefficient. Current numerical simulations give ass of order 10~^ or less, which 
would suggest that this coefficient is very small. On the other hand, these simulations have 
Va ~ Cs and are definitely not in the asymptotic regime where our scaling laws should 
be valid. We have already noted that the small difference in the exponent of Va in the 
dynamo growth rate and the dissipation rate, coupled to the presence of corrections to 
both these rates of order (Va/c^), makes it difficult to extrapolate from current results. If 
the saturation value of Va/cs approaches its asymptotic dependence on h/r gradually as 
/i/r — >■ then the final value of the coefficient will be much larger than 10~^. Bearing in 
mind the large role that numerical viscosity plays in the simulations ( [Vishniac 1995a| ) it 
seems prudent to regard the coefficient as an unknown numerical constant. 

On the other hand, since /i -C r for many realistic disks we can compare this dynamo 
mechanism to others based purely on the value of the exponent in the scaling relationship. 
Of course, a purely local mechanism co/ierent mechanism will not scale with (h/r) at all, 
although it might show some dependence on the local disk temperature. However, as we 
noted earlier this model seems to conflict with phenomenological studies of dwarf novae 
and x-ray transients. Internal waves, excited by tidal instabilities in binary system disks 
( [Goodman 1993D will produce an effective ass which scales as {h/rY ( [Vishniac, fc Diamond 



[1989[ ). This will be a competing mechanism for angular momentum transport in gas pressure 
dominated disks, and potentially the dominant one in radiation pressure dominated disks. 
(Although, such conditions are most likely in AGN disks, where the potential for the tidal 
excitation of waves is less certain.) Given the nonlocal nature of the angular momentum 
transport mediated by internal waves, the existence of a purely local mechanism might be 
important, even if it does not clearly dominate. When the disk is ionized and when internal 
waves are present, then the waves are capable of driving a dynamo with a growth rate 
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~ [h/rY^'^VL (IVishniac, Jin fc Diamond 1990 ) and perhaps faster, depending on the nature 
of the turbulent cascade of wave energy ( [Vishniac, fc Diamond 199^ ). The resulting value 
of ass will be ~ {h/rfl^{Pgas/P) ( |Vishniac, k Diamond 19'^^ , [Vishniac 1993BD . When 
these conditions are met this would appear to be a more important dynamo mechanism, 
although once again we note that nonlocal effects on the wave-driven dynamo make the two 
processes somewhat incommensurate. An equivalent estimate based on purely local physics 
was given by Meyer & Meyer- Hoffmeister (1983). However, this estimate is based on using 
large scale buoyant cells driven by magnetic buoyancy, a picture which is inconsistent with 
turbulence in the disk ([Zweibel fc Kulsrud 1975] , [Vishniac, fc Diamond 19921) . In addition, 
they assumed approximate isotropy of the helicity tensor and offered a calculation of a^r 
instead of age. This assumption of isotropy is inconsistent with the notion that the motions 
are driven by magnetic buoyancy, which for /i ^ r will have a time scale much longer than 
the local shearing time scale. Finally, we note that our result for the incoherent dynamo 
in an accretion disk is sensitive to our assumption that the process must be self-exciting. 
Given an external source of fluid turbulence, e.g. convection or the turbulent cascade 
of energy set off by finite-amplitude internal waves, the incoherent dynamo may have a 
significantly larger growth rate and give rise to a larger ass- 



4. Conclusions 



In this paper we have shown that mean-field dynamo theory allows for the existence 
of a new kind oi a — Vt dynamo, which we have named the incoherent a — Vl dynamo, in 
which there is no coherent helicity whatsoever. In this class of dynamos the magnetic field 
is driven by a combination of a random walk for and its shearing, which creates Bq. The 
resultant large scale field derives its organization from coherent shearing effects, rather than 
any loss of mirror symmetry in the turbulence. Although this kind of dynamo necessarily 
includes spontaneous field reversals, such reversals may occur at a rate which is some 
fraction of the dynamo growth rate. The existence of a mean-field dynamo in a fiow with 
a mean helicity of zero is interesting for its own sake, since it provides an example of how 
large scale order in the magnetic field can arise from the interaction between a large scale 
shear and statistically symmetric local motions. In this sense it represents an alternative 
to models which seek to explain dynamo activity through asymmetric turbulence and a 
coherent helicity. It differs from previous attempts to do without a coherent helicity (e.g. 
Montgomery et al. (1984) and Gilbert et al. (1988)) in that it does so without appealing 
to the other terms in equation ([^). This model is particular interesting in light of previous 
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claims (|Motfatt 1979 ) that the coherent a effect does not converge (although see Kraichnan 



(1979) for a counterargument). 

This dynamo is particularly interesting in light of simulations of magnetic field 
instabilities in accretion disks. We have suggested that this dynamo can operate successfully 
in accretion disks, but only to produce axisymmetric large scale fields. Comparing the 
growth rate for this dynamo with the buoyant loss rate for magnetic flux, we see that ^/ 
this is the only dynamo associated with magnetic shearing instabilities, then the large scale 
magnetic field will saturate when Va ~ {h/r)cs and ass ~ {H/rY. This result may seem 
somewhat odd, since the dynamic equations do not depend on r at all. However, the factor 
of r comes in through geometrical considerations, i.e. from considering the number of 
independent eddies in an axisymmetric magnetic field domain. In that sense it refers to the 
circumference of such an annulus rather than its radius. Consequently, when comparing 
numerical simulations to the predictions of this model one should substitute the azimuthal 
extent of the simulations for r. For current simulations this gives {h/r) ~ (27r)^^. In fact, 
since the turbulent eddies are longer in the azimuthal direction, the number of independent 
eddies that can be stacked end to end in current simulations is one or two, so that effectively 
{h/r) ~ 1. A critical test of this model would be to run simulations where this ratio was 
small and look for a drop in Va/cs in the saturated state. It's interesting to note that the 
only numerical simulation with no imposed field and with disk vertical structure does seem 
to have a coherent component to the helicity, which may have the wrong sign to drive a 
conventional dynamo. This may account for the large rate of spontaneous field reversals 
when vertical structure is put in the simulations. The existence of such a component is 
consistent with the existence of an incoherent a — VL dynamo, but only if its amplitude 
scales steeply with the strength of turbulence in the disk. 

Finally, we note that this model successfully reconciles phenomenological models of 
stellar accretion disks and the existence of a dynamo effect in a magnetized disk. The only 
drawback is that this model gives a relationship between ass and {h/r) which is probably 
too steep, implying the existence of other, more efficient dynamo mechanisms in accretion 
disks in binary systems. 
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Fig. 1. — Contours of the Be field in a space-time diagram for the one dimensional spatially 
and temporally incoherent dynamo model. 




Fig. 2. — A snapshot of the magnetic field and hehcity for the one dimensional spatially 
incoherent dynamo model. 




Fig. 3. — The evolution of a^g (normalized by the product of the rms velocity and {flTeddy) ^) 
in the upper half plane of run C from Brandenburg et al. (1995). 



